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FOREWORD 


This  document  Is  Volume  XVIII  of  the  Interim  Report  series  for  the  Passive 
Nosetip  Technology  (PANT)  prog ran.  A summary  of  the  documents  in  this  series  pre 
pared  to  date  is  as  fol’ows: 

Volume  1 - Program  Overview  (U) 

Vokv#s  H - environment  and  Material  Response  Procedures  for  Nosetip 
Design  (U) 


Volume  ?II 


Volume  IV 


Volume  V 


Volume  VI 
Volume  VII 

Volume  VIII 


Volume  IX 


- Surface  Roughness  Data 

Part  I - Experimental  Data 

Part  II  - Roughness  Augmented  Heating  Data  Correlation  and 
Analysis  (U) 

Part  III  - Boundary  Layer  Transition  Data  Correlation  and 
Analysis  (U) 

- Heat  Transfer  and  Pressure  Distributions  on  Ablated  Shapes 
Part  I Experimental  Data 

Part  II  - Data  Correlation 

- Definition  of  Shape  Change  Phenomenology  from  Low  Temperature 
Ablator  Experiments 

Part  I - Experimental  Data,  Series  C (Preliminary  Test 
Series) 

Part  II  -Experimental  Data,  Series  D (Final  Test  Series) 

- Graphite  Ablation  Data  Correlation  and  Analysis  (U) 

- Computer  User's  Manual,  Steady-State  Analysis  of  Ablating 
Hosetlps  (SAA.TT)  Program 

- Computer  User's  Manual,  Passive  Graphite  Ablating  Nosetip 
(PAGAN)  Program 

- Unsteady  flow  ol  Ablated  Nosetip  Shapes  - PANT  Series  G Test 
and  Analysis  Repurt 
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Volume  X - Summary  of  Experimental  and  Analytical  Results  for  the  Period 
Nay  1973  to  OecMter  1974. 


Volts*  XI  - Analysis  and  Revlow  of  the  ABRES  Combustion  Test  Facility 

for  High  Pressure  Hyperthermal  Reentry  Nosetlp  Systems  Tests 

Volume  XII  - Nosetlp  Transition  and  Shape  Change  Tests  In  the  AFFDL  50  MM 
RENT  Arc  - Data  Report 

Volume  XIII  - An  Experimental  Study  to  Evaluate  Heat  Transfer  Rates  to  Scal- 
loped Surfaces  - Data  Report 


Volume  a.V 


Volume  XV 
Volume  XVI 
Volume  XVII 


- An  Experimental  Study  to  Evaluate  the  Irregular  Nosetip  Shape 
Rig i me  - Data  Report 

- Roughness  induced  Transition  Experiments  - Data  Report 

- Investigation  of  Erosion  Mechanics  on  Reentry  Materials  (U) 

- Computer  User’s  Manual.  Erosion  Shape  (EROS)  Computer  Program 
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ABSTRACT 


The  results  of  noset  ip  response  calculations  using  the  newly-developel,  EROsior 
Shape  (EROS)  computer  -ode  are  docuaented,  The  code  er^oles  the  evaluatior  of  nosetio 
coupled  erosion  and  ablation  and  shape  change  for  arbitrary  materials  including  fine 
weave  carbon/carbon,  polycrystalline  graphite,  and  reinforced  phenolics.  Transient 
heat  conduction  Is  also  treated  through  an  Innovative  explicit  differencing  scheme. 
Calculations  are  presented  for  wi» d tunnel  test  environments,  ICSK  ciear  air  flight  con- 
ditions, and  Terrier-Recruit  sounding  rocket  flight  conditions.  In  the  analysis  matrix, 
modeling  techniques  are  varied  to  Identify  response  sensitivity  to  modeling  uncertain- 
ties. 

The  primary  conclusions  are  that: 

e Transition  altitude  and  overall  recession  of  ICBM  nosetips  in  clear  air  are 
critically  sensitive  to  meterlal  roughness  characteristics. 

• Hydrometeor  Impacts  cause  significant  nosetip  erosion  ■i»ass  loss  and  increased 
surface  thermochemical  mass  loss. 
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INTRODUCTION 

Nosetip  ablation  and  erosion  performance  is  critical  to  the  survival  of  high 
performance  1CBM  reentry  veht-le  systems.  Analytical  procedures  have  been  developed 
under  the  Passive  Nosetip  Technology  t PANT ) program  (Contract  E04/01 -71-C-OOP7 ) and 
the  Hlnuteman  Natural  Hazards  Assessment  Program  (Contract  T 0470 1 -74-C-0069)  to 
evaluate  nosetip  design  concepts  and  to  assess  survival  probabilities.  The  most 
recent  nosetip  shape  and  thermal  analysis  code  includes  the  hydrometeor  effects 
modeling  described  in  Reference  1 and  a transient  heat  condu' tion  package  described 
in  Reference  2.  The  code  is  entitled  EROsion  Shapo  (EROS)  computer  program,  and 
Reference  3 is  the  computer  user's  manual.  Results  generated  during  the  prelimi- 
nary exercise  of  the  code  are  presented  herein 

Solutions  have  been  obtained  for  conditions  corresponding  to: 

• Pressure,  calorimeter,  and  shape  change  tests  from  PANT  wind 
tunnel  series 

a Sounding  rocket  erns mn-aoiat ion  flight  environments 

a ICBM  flight  environments 

The  comparisons,  included  herein,  between  data  and  code  results  for  a wide  range  of 
modeling  assumptions  and  variations  provide  insight  into  the  validity  and  consis- 
tency of  the  computational  tool.  The  majority  of  solutions  were  performed  for  a 
material  that  is  well  characterized  and  typical  of  advanced  nosetip  design.  Shape 
change  solutions  were  also  performed  for  a SAMS  carbon/carbon  composite  nosetip  and 
a SAMS  carbon  phenolic  nosetip. 

Section  2 presents  a brief  summary  of  me  EROS  code  features.  Clear  air 
flight  and  ground  test  solutions  showing  response  sensitivity  to  modeling  pertur- 
bations are  given  in  Section  3,  and  weather  related  sensitivity  analvses  are  presented 
in  section  4.  Carbon/carbon  composite  and  carbon  pherolic  weather  response  solu- 
tions are  presented  in  Section  5,  and  conclusions  are  summarized  in  Section  6. 
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SfCTION  2 


ABLATION/SHAPE  CHANGE  ANAtYSIS  PROCEDURES 

The  EROsion  Shape  (EROS)  computer  code  confines  the  environmental  modeling 
techniques  developed  under  the  PANT  and  Natural  Hazards  Assessment  Program  with 
the  in-depth  transient  conduction  routines  developed  at  the  Aerospace  Corporation. 
This  code  numerically  models  the  shape  history  of  an  ablating  nosetip  being  sub- 
jected to  a reentry  environment.  The  code  calculates  the  inviscid  flow  and  heat 
transfer  distribution  for  arbitrarily-shaped,  axisymmetric,  blunt  bodies  in  hyper- 
sonic flow.  In  addition,  the  boundary  layer  and  heat  transfer  distributions  are 
modeled  for  a variety  of  environments  including  the  effects  of  hydrometer  erosion. 
Two  material  response  prediction  techniques  are  available: 

• Steady  state  surface  energy  balance 

• A surface  energy  balance  fully  coupled  to  the  in-depth 
transient  thermal  response 

The  in-depth  thermal  response  is  capable  of  calculating  the  three-dimensional 
temperature  field  and  surface  recession  of  nosetips  at  angle  of  attack.  However, 
since  the  environment  package  is  limited  to  axisymmetric  geometries,  the  present 
code  is  restricted  to  nosetips  at  zero  angle  of  attack.  Three  additional  assump- 
tions implicit  in  the  formulation  of  these  computerized  ablation/shape  change 
analyses  procedures  are  as  follows: 

• Flow  separations  do  not  occur 

• Oscillatory  inviscid  flow  fields  do  not  develop 

• Secondary  shocks  do  not  occur 

The  elements  of  shape  change  analyses  in  the  EROS  code  are  illustrated  in  the 
computation  procedure  sketch  in  Figure  2-1.  Briel  descriptions  of  each  of  these 
elements  are  presented  in  tl  s section.  The  inviscid  flow  field  modeling  is 
described  in  Section  2.1.  Boundary  layer  analysis  procedures  are  discussed  in 
Section  2.2.  The  graphite  ablation  model  i.;  sunmarized  in  Section  2.3  with  a 
description  of  the  body  movement  and  shape  redefinition  procedures  in  Section  2.4. 
Erosion  modeling  capabilities  are  reviewed  in  Section  2.5. 
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2.1 


INVISCID  FLOW  FIELD 


Calculation  of  the  inviscid  flow  field  over  an  ablating  shape  provides  the 
required  boundary  conditions  for  the  boundary  layer  solution.  One  aspect  of  the 
inviscid  flow  field  solution  in  the  PANT  computer  codes  is  the  pressure  distribu- 
tion calculation.  The  pressure  distribution  over  a nosetip  contour  is  calculated 
based  on  regional  correlations,  defined  as  (1)  subsonic  forebody,  (2)  supersonic 
forebody  including  recompression,  and  (3)  the  aft  body.  The  pressure  distribution 
in  the  subsonic  region  depends  upon  the  free  stream  Mach  number.  The  correlation 
in  this  region  is  an  empirical  extension  and  synthesis  of  the  modified  Newtonian 
correlation,  valid  for  spheres,  but  Including  a correlation  for  flat  faced  cylin- 
ders. This  correlation  strongly  depends  on  the  sonic  pyint  location  which  is 
evaluated  assuming  modified  Newtonian  flow  with  perturbations  in  the  nominal  solu- 
tion depending  on  (1)  free  stream  Mach  number,  (2)  specific  heat  ratio,  (3)  nose- 
tip  bluntness,  and  (4)  surface  streamline  recompression  on  biconic  shapes. 

On  the  supersonic  forebody  of  the  nosetip,  pressure  distributions  are  com- 
puted either  using  the  modified  Newtonian  expression,  or  for  biconic  type  confi- 
gurations, using  a conic  surface  recompression  correlation.  The  cone  recompres- 
sion model  is  based  on  sphere/cone  and  ellipsoid/cone  exact  solutions  perfort oo 
at  various  Mach  numbers. 

The  correlation  for  aft  cone  pressures  is  one  developed  at  the  Aerospace 
Corporation  (Reference  4).  Its  form  is  such  that  the  predicted  conic  pressure 
asymptotically  approaches  the  sharp  cone  value.  Coupling  the  three  regimes  into 
a continuous  distribution  is  achieved  by  a smoothing  technique  based  on  a weighted 
average  of  the  incremental  modified  Newtonian  expression  and  a linear  decay 
expression.  Details  regarding  the  various  aspects  of  the  pressure  distribution 
correlations  are  available  in  Reference  3. 

2.2  BOUNDARY  LAYER  NEAT  AND  MASS  TRANSFER 

The  boundary  layer  heat  and  mass  transport  events  are  mooeled  using  a film 
coefficient  approach.  The  two  aspects  of  the  heating  distribution  predictions 
which  were  emphasised  in  the  upgrading  of  the  EROS  code  were  (1)  transition  and 
transitional  flow  and  (2)  rough  wall  heatirn  perturbations  to  both  laminar  and 
turbulent  flow.  Development  of  the  transition  model  is  discussed  in  Reference  5. 
The  laminar  and  turbulent  rough  wall  heating  augmentation  models  are  discussed  in 
Reference  6.  The  sensitivity  of  ablation  shape  change  predictions  for  ATJ-S 
graphite  to  uncertainties  in  the  rough  wall  heating  levels  tna  transition  location 
are  reviewed  in  Reference  7. 
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2.3  THERMOCHEMICAL  ABLATION  MODEL 


In  the  EROS  computer  code,  the  ablation  thermochemistry  events  are  modeled 
by  the  generalized  approach  described  in  Reference  8.  Through  a tabular  coupling 
to  a thtrmc  hemistry  computer  cede,  the  surface  energy  balance  routine  can  consider 
any  combination  of  surface  materiel  elemental  commltions  in  an  air  environment. 

The  film  coefficient  approach  enables  the  modeling  of  heterogeneous  reaction  and 
sublimation  kinetics,  unequal  species  diffusion  coefficients,  and/or  unequal  heat 
and  mass  transfer  coefficients.  For  the  case  of  graphite  in  air,  dominant  surface 
ablation  and  reaction  products  Include  CO,  CNr  C?N,  C2N2,  NO,  N02,  C02,  and  carbon 
vapor  species. 

2.4  BODY  MOVEMENT  Niu  SURFACE  SMOOTHING 

A significant  amount  of  development  work  in  the  area  of  shape  stability  and 
body  point  movement  has  occurred  within  the  PANT  program.  Previous  stability  problems 
resulted  from  surface  movement  schemes  which  allowed  body  point  migration  toward  the 
stagnation  point.  The  steady  state  option  of  the  EROS  code  eliminates  these  problems 
by  restricting  body  point  movement  to  be  along  lines  of  constant  radius,  in  addition, 
an  effective  nose  radius  scheme  is  used  in  the  shape  change  procedure.  This  specifi- 
cation controls  only  the  detailed  shape  dependent  environmental  parameters  at  the 
stagnation  point  and  does  not  limit  or  redefine  the  overall  nosetip  shape.  These 
shape  change  mechanics  and  body  point  movement  criteria  are  fundamental  to  the  com- 
putational stability  of  the  EROS  code. 

2.5  EROSION  MODELING 

Coupled  erosion  ablation  computations  '.wy  be  performed  with  the  EROS  code. 
Surface  erosion  due  to  hydrometeor  impacts  is  currently  modeled  for  graphitic  type 
brittle  materials  and  malleable  type  metal  materials.  Various  erosion  modeling 
assumptions  available  in  the  EROS  code  are  docunented  in  Reference  1 and  summarized 
in  Section  4. 
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SECTION  3 


ABLATION  AND  SHAPE  CHANGE  MODELING  SENSITIVITY  STUDIES 

The  EROS  code  w<>s  exercised  to  determine  code  reliability  and  computational 
consistency  for  noneroding  environments  including  wind  tunnel  tests  and  clear  air 
flight  environments.  Solutions  were  generated  for  reasonable  perturbations  to 
several  modeling  features  which  have  recently  been  upgraded  and  incorporated  in 
the  PANT  codes  (in  this  case,  the  EROS  code).  From  these  solutions,  the  sensiti- 
vity of  the  EROS  code  predictions  to  key  modeling  uncertainties  was  demonstrated, 
and  code  numerical  difficulties  were  identified. 

In  Section  3.1  comparisons  with  PANT,  wind-tunnel  test  data  are  presented. 
Shock  shape  and  heat  transfer  distributions  over  ablated-shape  calorimeters  are 
shown,  and  the  results  cf  two  low- temperature-ablator  shape  change  predictions 
are  cor, pared  to  data.  In  Section  3.2,  several  code  calculations  for  two  typical 
reentry  cases  are  presented.  For  each  flight  case,  a baseline  steady  state  abla- 
tion solution  is  presented.  The  effects  of  perturbed  modeling  are  then  referenced 
to  the  baseline  result. 

3.1  COMPARISONS  TO  PANT  WIND  TUNNEL  TEST  DATA 

Wind  tunnel  test  data  from  calorimeter  and  low-temperature-ablator  shape 
change  models  provide  data  to  assess  the  nosetip  response  methodology  in  the  EROS 
code.  For  example,  ablated  shape  calorimeter  data  from  the  Senes  B tests 
(Reference  9)  and  the  more  recent  Series  J tests  (Reference  10)  are  ideal  fcr 
evaluating  shock  shape  and  heat  transfer  modeling  techniques.  In  addition,  exact, 
inviscid  flow  field  solutions  have  been  performed  for  the  Series  B test  geometries 
at  the  Mach  5 wind  tunnel  condition  and  at  a Mach  20  condition.  The  high  Mach 
number  exact  solutions  enable  comparison  at  flight  conditions  where  data  are 
unavailable. 

Four  model  configurations  were  selected  for  analysis  of  EROS  shock  shape 
and  heat  transfer  modeling.  These  have  nosetip  profiles  characterized  as: 

• Convex  (Series  B) 

• Concave  (5eries  B) 
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• 45°  conic  (Series  B) 

• 60°  conic  (Series  J) 

In  each  case,  a 0. 4-inch  radius  stagnation  region  fairs  into  the  nosetlp  profile 
which  fairs  into  an  8°  aft  cone.  The  model  profiles  are  shown  along  with  the 
computed  and  measured  shock  shapes  in  Figures  3-1  through  3-4. 

T!  e EROS  modeling  of  shock  shape  (Reference  1)  includes  two  improvements 
over  the  previous  PANT  model  (Reference  11).  First,  an  alternate  shock  standoff 
correlation  is  used  and,  second,  an  improved  procedure  is  used  for  computing  the 
shock  shape  downstream  o*  the  sonic  flow  line.  No  change  was  made  to  the  sub- 
sonic region  shock  shape  computation  procedure.  The  convex  shape  comparisons  in 
Figure  3-1  demonstrate  a large  Improvement  in  the  shock  standoff  calculation  and 
a small  improvement  in  overall  shape.  The  comparisons  for  the  concave  pit>file  in 
Figure  3-2  show  a substantial  improvement  in  the  shock  sh4>e  both  at  Nach  5 and 
20,  whereas  Figure  3-3  shows  no  difference  in  the  results  for  the  45°  biconic. 

Both  the  previous  model  and  the  EROS  model  produced  an  accurate  shock  shape  for 
the  45°  conic.  The  results  for  the  60°  conic  profile  in  Figure  3-4  are  not  as 
good,  however.  The  centerline  shock  standoff  calculation  is  closer  to  the  data, 
but,  since  the  subsonic  region  model  was  not  changed,  the  overall  shock  shape  in 
front  of  this  high  angle  conic  was  not  changed  from  the  previous  result.  No  high 
Mach  number  solution  was  run  for  the  60°  conic. 

Smooth  wall  heat  transfer  distribution  predictions  for  the  convex,  concave, 
and  45  degree  conic  configurations  are  compared  to  SAANT  code  predictions  (Refer- 
ence 12)  and  data  in  Figure  3-5.  Of  the  three,  only  the  45°  conic  prediction 
differs  substantially  from  the  SAANT  results.  The  EROS  code  incorporates  the 
composite  heat  transfer  prediction  technique  described  in  Reference  3 and  in- 
cludes different  modelings  of  boundary  layer  growth,  entropy  layer  swallowing, 

Reynolds  analogy,  and,  of  course  rock  shape.  For  the  convex  and  concave  shape, 
the  modifications  compensated  to  give  sinilar  predictions. 

As  an  exercise  of  the  shape  change  numerical  procedures  in  the  EROS  com- 
puter code,  two  camphor  shape  change  solutions  were  generated  for  comparison  with  PANT 
Series  D Runs  207  and  208.  The  model  in  Run  207  was  tested  at  a relatively  high  Rey- 
nolds number  condition  and  experienced  shape  change  to  a blunt  turbulent  configuration. 
Run  208,  performed  at  a lower  Reynolds  number  condition,  resulted  in  a slender 
shape.  Figures  3-6  and  3-7  present  the  shape  and  recession  predictions  for  Runs 
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Figure  3-2.  Shock  shape  comparisons  for  concave  nosetip  profile. 
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Figure  3-5.  Smootfcwall  heat  transfer  distribution  comparisons. 


(b)  Stagnation  point  recession. 

Figure  3-6.  Run  207  camphor  shape  change  prediction 
( Rea,  « 10  x 10‘/ft). 
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207  and  208,  respectively.  When  compared  with  the  wind  tunnel  data,  one  observes 
that  the  predicted  shape  generally  agrees  with  data.  However,  the  solution  for 
Run  208  underpredicts  the  stagnation  point  recession.  This  undcrprediction  is  at- 
tributed to  the  uncertainty  in  forecone  macroroughness  modeling. 

3.2  REENTRY  VEHICLE  NOSETIP  SOLUTIONS 

The  upgraded  PANT  computer  code  (EROS)  was  exercised  for  two  high  perfor- 
mance  reentry  situations.  The  objectives  of  the  calculations  were  two  fold: 

• Exercise  and  "debug"  erosion/ablation  coupling  logic  for  clear  air 
solutions. 

• Understand  the  sensitivity  of  clear  air  nosetip  shape  and  recession 
predictions  to  uncertainties  in  and  modifications  to  modeling 
techniques. 

Trajectory  parameters  and  nosetip  configurations  for  these  two  cases  are  described 
in  Section  3.2.1.  In  Section  3.2.2,  the  matrix  of  clear  air  response  modeling 
variations  are  summarized.  Results  are  discussed  in  Section  3.2,3. 

3.2.1  Trajectory  Parameters  and  Nosetip  Configurations 

Two  typical  trajectory  and  nosetip  configuration  cases  were  analyzed. 

Table  3-1  sumnarizes  basic  features  of  the  cases.  Trajectory  properties  (vehicle 
velocity  and  stagnation  pressure)  are  plotted  versus  attitude  and  reentry  time  for 
the  two  cases  in  Figure  3-8.  Initial  configurations  are  shown  in  Figure  3-9. 


TABLE  3-1.  TRAJECTORY/CONFIGURATION  PARAMETERS 


Case  1 

Case  2 

Entry  Velocity,  kft/sec 

22.8 

22.3 

Entry  Angle,  Degree 

-38 

-26 

Ballistic  Coefficient,  lb^/ft2 

2150 

3300 

Nose  Radius,  in. 

1.25 

0,75 

Cone  Half  Angle,  Degree 

6.7 

6.3 

Material  Type 

Fine  Grain  Graphite 

Polycrystalline  Graphite 
or  Tine  Weave  Carbon/Carbon 

Nosetip  Design 

Shell 

Plug 
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b)  Case  2 (R^  = 0.75  in.) 

Figure  3-9.  Initial  configurations  for  clear  air  flight 
response  sensitivity  calculations. 
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The  two  cases  represent  a significant  variation  in  several  parameters, 
yet  within  current  conditions  of  interest  to  nosetip  and  reentry  vehicle  system 
designers.  Both  nosetip  materials  are  basically  graphitic  with  the  primary  dis- 
tinction (from  a shape  change  point  of  view)  coming  from  the  micro  or  intrinsic 
roughness  used  as  the  baseline  for  analysis  sensitivity  studies,  The  smaller 
nose  radius,  high  ballistic  coefficient  case  (Case  2)  is  representative  of  a 
higher  recession  situation  than  Case  1.  Thus  the  two  cases  exercised  the  shape 
code  differently. 

3.2.2  Matrix  of  Ablation  Modeling  Parametric  Solutions 

For  each  trajectory-configuration  case,  baseline  modeling  and  material 
property  values  were  selected  and  a solution  was  generated.  Systematic  variations 
from  the  baseline  were  then  defined  and  the  cases  were  run  for  comparison. 
Variations  in  the  following  were  considered: 

• transition  due  to  microroughness  variations 

§ transition  due  to  alternate  transition  criteria  specifications 

• turbulent  heat  transfer  due  to  macroroughness  (scallop)  effects 

• shock  shape  definition  procedures 

Table  3-2  summarizes  the  modeling  assumptions  used  for  the  respective  trajectory- 
configuration  calculations.  Shape  profile  histories  for  the  fifteen  solutions 
are  shown  in  the  Appendix.  Suitable  comparisons  between  respective  solutions  are 
presented  in  the  following  section. 

3.2.3  Results  of  Parametric  Solutions 

Parametric  calculations  including  variations  in  critical  modeling  provide 
an  excellent  basis  to  check  out  computer  code  logic  and  simultaneously  to  clarify 
the  effects  of  modeling  uncertainties  on  computations.  In  this  section,  several 
results  from  the  parametric  analysis  matrix  (Table  3-2)  are  compared  and 
conclusions  are  summarized.  For  reference,  the  nominal  baseline  solutions  for 
cases  1 and  2 are  shown  in  Figure  3-10  and  3-11,  respectively.  The  figures 
include  both  the  shape  profile  history  (shape  profile  at  roughly  10  kft  inter- 
vals) and  the  stagnation  point  recession  history  for  the  two  cases.  Both  the 
shape  and  recession  responses  are  significantly  different. 

Two  solutions,  in  addition  to  the  baseline  runs,  were  performed  for  both 
cases  1 and  2 to  assess  the  affect  of  material  microroughness  on  the  nosetip 
transition  onset  altitude  and  on  the  resulting  stagnation  point  recession. 
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TABLE  3-2.  MATRIX  FOR  ABLATION  MODELING  PARAMETRIC  SOLUTIONS 
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Figure  j-10.  Nominal  baseline  prediction 
{Solution  1)  for  Case  1. 
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(b)  Stagnation  point  recession. 

-11.  Nominal  baseline  prediction  (Solution  6) 
for  Case  2. 
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Recession  histories  for  these  six  predictions  (Solutions  1,  2,  3,  6,  7,  and  8 In  Table 
3-2)  are  compared  In  Figure  3-12.  For  each  case,  larger  roughness  means  earlier 
transition  and  more  total  recession.  Some  amount  of  time  Is  required  between 
transition  onset  and  the  time  when  stagnation  point  recession  Is  affected.  During 
this  period,  the  nosetip  shape  Is  changing  from  the  blunt,  laminar  profile  to  a 
conic  type  shape.  Once  the  sharpening  process  nears  completion,  the  stagnation 
region  radius  of  cur\'tura  decreases,  and  heating  and  recession  rate  increase.  The 
earlier  transition  onset  occurs,  then  the  sooner  the  shape  will  become  sharp.  Figure 
3-1 2b  indicates  the  significant  Increase  In  recession  rate  once  sharpening  Is 
complete. 

Figures  3-13  and  3-14  further  demonstrate  the  sensitivity  of  the  clear  air 
nosetip  response  prediction  to  the  value  of  the  material  microroughness.  Figure 
3-13  gives  nosetip  transition  onset  altitude  as  a function  of  microroughncss. 

Symbols  indicate  the  results  from  the  parametric  solution  matrix,  including  those 
solutions  using  the  updated  PANT  transition  criterion.  The  figure  indicates  the 
following: 

• Transition  onset  altitude  is  very  sensitive  to  microroughness; 

• Essentially  no  difference  exists  between  cases  1 and  2 (effects  of 
nose  radius  and  trajectory  differences  compensate); 

• Use  of  the  updated  PANT  transition  criterion  gives  a slightly  higher 
transition  onset. 

Figure  3-14  shows  the  sensitivity  of  overall  stagnation  point  recession  to  the 
material  microroughness  value  for  cases  1 and  2.  The  estimated  trends  art*  based 
on  the  anticipated  sharpening  altitude  and  recession  rate.  An  additional  data 
point  from  Reference  13  has  been  added.  Note  the  significant  increase  in  reces- 
sion that  can  occur  for  a small  change  in  microroughness. 

Once  transition  onset  occurs,  the  turbulent  boundary  layer  interacts  with 
the  ablation  process  to  increase  the  surface  roughness  through  the  formation  of 
crosshatches  or  scallops  on  the  nosetip.  Based  on  results  from  Reference  14,  the 
effective  scallop  roughness  (macroroughness)  has  been  modeled  as  a function  of 
edge  pressure, 

Kt  = APe'°v77  t3*1 ) 

where  A has  been  defined  for  materials  of  interest  by  evaluating  ablation  data  taken 
in  the  AFFDl.  50  MW  arc  jet.  A convenient  means  of  varying  the  macroroughness  within 
the  uncertainty  of  the  50  MW  data  analysis  is  to  use  the  stagnation  pressure  in  place 
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of  the  edge  pressure  In  the  above  correlation.  This  was  done  for  solutions  4 and 
11  of  the  parametric  matrix  (Table  3-2)  as  shewn  in  Figure  3-15,  this  modeling 
change  has  the  following  effects: 

• Sharpening  time  is  longer  because  the  turbulent  forecone  heating 
and  recession  rate  are  less; 

• Recession  rate  after  sharpening  is  roughly  12  percent  lower  when 
stagnation  pressure  is  used  (See  Case  2,  Figure  3-1 5b ) . 

Nosetip  response  predictions  comparing  the  SAANT  shock  model  (Reference  11) 
and  the  revised  shock  model  (Reference  1^  demonstrated  the  effect  of  this  modeling 
update.  Figure  3-16  shows  the  final  shape  profile  comparison  for  Case  1 and  the 
recession  history  comparison  for  Case  2.  For  the  large  nose  radius  case,  which 
experiences  relatively  little  turbulent  recession,  the  effect  was  negligible. 

For  the  small  nose  radius  case,  the  new  shock  model  produced  slightly  more  stagna- 
tion point  recession  and  a slightly  highe*”  angle  turbulent-forecone  shape.  This 
trend  agrees  with  the  expected  effects  of  the  improved  shock  model  since  the 
principle  modification  came  in  the  supersonic,  ablated  shape  situations,  not  in 
the  laminar  or  sharpening  situations. 

The  importance  of  modeling  transitional  heating  was  assessed  by  a compari- 
son between  the  baseline  solution  and  a solution  based  on  a "sharp"  transition 
assumption  (Solutions  6 and  10  in  Table  3-2).  Sharp  transition  means  immediate 
increase  to  fully  turbulent  heating  rate  at  the  nosetip  location  where  transition 
is  predicted.  Figure  3-17  shows  the  recession  history  comparisons  for  Case  2 
(Rn  = 0.75  in.).  The  primary  effect  of  this  modeling  variation  was  a reduction 
in  the  time  required  for  the  nosetip  to  sharpen  from  the  shape  at  transition 
onset  to  a fully  turbulent  shape  where  recession  increases  significantly.  For 
the  small  nose  radius  analysis  case,  the  net  recession  increased  roughly  10 
percent  due  to  the  higher  sharpening  altitude. 

Two  nosetip  response  predictions  for  the  small  ncse  radius  case  were 
obtained  using  the  transient,  2-0  heat  conduction  logic  in  the  EROS  code.  Com- 
parable steady  state  solutions  indicated  that  modeling  of  transient,  2-D  heat 
conduction  was  important  but  not  critical  to  the  evaluation  of  shape  and  recession 
trends.  Two-dimensional  heat  conduction  towards  the  nosetip  centerline  generally 
increases  stagnation  point  recession  and  produces  a slightly  larger  nose  radius 
shape  than  predicted  in  steady  state.  A secondary  effect  of  the  stagnation  region 
blunting  can  be  a reduction  in  heat  transfer  and  a net  reduction  in  recession. 
These  trends  can  be  seen  in  the  transient,  steady  state  comparisons  in  Figure  3-18, 
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(b)  Stagnation  point  recession  (Case  2). 


Figure  3-16.  Effect  of  shock  model  upon  nosetip  response 
prediction. 
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Figure  3-17.  Effect  of  transitional  heating 

model  on  Case  2 recession  history. 
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Figure  3-18.  Effect  of  transient  heat  conduction 
on  Case  2 recession  history. 
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In  summary,  several  aspects  of  modeling  were  addressed  in  the  parametric 
code  check  out  matrix.  The  solutions  confirmed  that  the  code  was  functioning 
properly  in  the  reconmended  mode  of  operation.  The  solutions  indicated  that  the 
most  important  code  input  is  the  material  microroughness  since  this  parameter 
controls  transition  onset  attitude.  The  response  predictions  were  generally  less 
sensitive  to  other  modeling  changes  and  uncertainties. 
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SECTION  4 

WEATHER  EFFECTS  SENSITIVITY  SOLUTIONS 

The  EROsion  shape  (EROS)  computer  code  combines  the  environmental  and  shape 
modeling  techniques  developed  under  the  PANT  program  with  weather  effects  modeling 
developed  under  the  basic  National  Hazards  Assessment  Program  and  the  preliminary 
studies  under  this  effort  (Reference  1).  In  this  section,  the  results  of  the 
initial  code  solutions  generated  using  the  weather  related  modeling  features  are 
presented.  These  solutions  were  performed  primarily  to  check  out  code  logic. 

However,  the  resulting  sensitivity  study  does  provide  a basis  for  identifying 
critical  modeling  assumptions. 

Two  sounding  rocket  test  flights  performed  under  the  SAMSO  sponsored  SAMS/ 
Terrier-Recruit  series  were  selected  for  analysis.  These  relatively  low  velocity 
flights  of  graphitic  type  nosetips  provided  weather,  trajectory,  and  nosetip  response 
data  for  code  assessment  purposes.  Trajectory,  configuration,  and  weather  parameters 
for  the  two  flights  are  presented  in  Section  4.1, 

Each  of  the  weather  effects  modeling  features  developed  in  this  program  and 
reported  in  Reference  1 was  exercised  in  the  check  out  matrix.  These  included: 

• nominal  and  lower  bound  mass  loss  correlations  for  the  flight  material 

« crater  formation  and/or  crater  healing  effects  on  rough  wall  heat 

transfer 


• erosion  related  heating  augmentation 

• particle  slowdown  and  mass  loss  in  the  shock  layer. 

• particle  stripping  and  breakup  in  the  shock  layer 

The  specific  analysis  matrix  is  presented  in  Section  4,2. 

Comparisons  of  results  to  each  other  and  to  the  nosetip  response  data  are 
presented  in  Section  4.3. 


4.1  FLIGHT  TEST  CONDITIONS 

Two  similar  flight  tests  were  analyzed  using  the  ERCS  computer  code.  The 
tests  were  conducted  from  Wallops  Station,  Wallops  Island,  Virginia  in  February 
and  March  of  1972.  The  flights  were  performed  by  Sandia  Corporation  under  joint 
sponsorship  of  Sandia  and  SAMSO.  The  booster  for  each  flight  was  the  SAMS/Terrier 
Recruit  which  achieved  peak  velocities  of  roughly  8600  ft/sec.  A recovery  package 
enabled  retrieval  of  the  nosetip  payload.  The  flights  are  designated  Sandia  Test 
No.  R341412  and  R341413  or,  more  typically,  SAMS  6 and  SAMS  7.  Trajectory,  con- 
figuration, and  nosetip  response  -lata  are  given  in  Reference  15.  Weather  data 
were  obtained  from  Reference  16.  Key  flight  test  parameters  are  sumnai  ized  in 
Table  4-1.  Pertinent  trajectory  quantities  are  shown  in  Figure  4-1;  the  initial 
configuration  for  both  tests  in  shown  ir,  Figure  4-2,  and  the  weather  parameters 
(liquid  water  content  and  median  particle  diameter)  used  in  the  computer  analyses 
are  shown  in  Figure  4-3. 

TABLE  4-1.  SAMS  FLIGHT  TEST  PARAMETERS 


SAMS  6 

SAMS  7 

Peak  Velocity,  ft/sec 

8602 

8575 

Peak  Stag,  Press.,  atm. 

50 

57 

Nose  Radius,  in. 

0.625 

0.625 

Cone  Angle,  degrees 

9.0 

9.0 

WSI,  kmJ-gm/m5 

7.5 

2.9 

4.2  EROSION  MODELING  SENSITIVITY  MATRIX 

For  each  flight  case  the  matrix  of  EROS  code  solutions  shown  in  Table  4-2 
were  generated.  The  various  models  are  discussed  in  detail  in  Reference  1 and 
summarized  below 

Hass  Loss  Correlations 

The  nominal  mass  loss  correlation  used  for  the  EROS  code  solutions  was  de- 
rived from  single  impact  and  ballistic  range  erosion  data.  The  form  of  the  cor- 
relation is 

G <*  V1  *a  sin,*‘28 
nom 
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TABLE  4-2.  EROSION  MODELING  SENSITIVITY  MATRIX 


where 


G„„.  = ratio  of  ejected  mass  to  Incident  mass 
non 

V * Incident  particle  velocity  at  the  surface 
0 * Incident  impact  angle 

This  "G-Law"  relation  Is  given  as  Equation  3-9  of  Reference  1. 

A second  mass  loss  correlation  used  to  assess  the  sensitivity  to  mass  loss 
was  defined  as  the  lower  bound  of  the  erosion  data.  The  ratios  of  the  lower  mass 
loss  correlation  to  the  nominal  Is: 


These  two  correlations  give  the  same  mass  loss  at  V = 1734  ft/sec  but  the  lower 
correlation  gives  roughly  30  percent  less  mass  loss  at  the  Terrier-Recruit  peak 
velocity  condition  (8600  ft/sec). 

Crater  Roughness  Heating 

The  modeling  of  crater  effects  on  surface  heat  transfer  includes  two  parts: 

e evaluation  of  equivalent  peak-to- valley  sandgrain  roughness  dimension 

• application  of  the  PANT  roughwall  heat  transfer  modeling. 

The  application  of  the  PANT  roughwall  correlation  is  described  in  Reference  11. 

The  crater  roughness  is  evaluated  from  the  erosion  mass  loss  and  an  assumed  crater 
configuration  assumption.  For  a hemispherical  crater,  the  crater  roughness,  rc, 
(Equation  3-16  of  Reference  1)  is: 


where 

G is  the  mass  loss  ratio 

Pp  is  the  particle  density 

Pt  is  the  surface  material  density 

d„  is  the  incident  particle  diameter 
P 
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The  solutions  which  Include  the  crater  effects  used  the  larg  of  either  the  crater 
roughness,  the  intrinsic  microroughness  In  laminar  flow  (0.4  mils),  or  the  scallop 
marroroughness  Is  turbulent  flow.  Scallop  dimensions  were  found  front 

k - 0.93  Pe-0*77 


For  one  matrix  point,  the  effects  of  crater  healing  were  modeled.  Healing  refers 
tc  the  reduction  of  crater  roughness  between  Impacts  due  to  ablation  from  the  top 
nr  the  crater.  The  healing  model  is  formulated  such  that  the  average  roughness 
dimension  between  successive  Impacts  Is  used  for  heat  transfer  calculations,,  lor 
example  if  the  time  between  successive  Impacts  (obscuration  time,  t ^s)  equals  the 
time  to  ablate  away  the  crater  dimension  (t^),  Then  the  average  crater  roughness 
is: 


k - rc/2 


tobs  “ tdes 


A1sc  for  W < tdes 


k * rc  - ~jr 


obs 


For  ‘obs  5 *des 


k ■ rca'2S‘obs 


where  s is  see 


tobs  < tdes 


t , > t . 

obs  des 

ion  rate  from  the  top  of  the  crater. 


Erosion  Augmentation 

In  the  erosive  environment,  the  presence  of  hydrometeors  and  surface  erosion 
ejecta  augment  the  heat  tr.  r to  the  nosetip.  The  physical  mechanisms  for  this 
heating  increase  are  not  . . defined,  but  the  effect  has  been  correlated  by  the 
following  expression  (Reference  17): 

* ■ <H0  - H„>  CH,0 


where  q is  incident  laminar  heat  flux 

HQ  - Hw  is  the  total  to  wall  enthalpy  potential 

p®u®  is  the  freestream  air  mass  flux 
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c 


H,u 


Is  the  erosion  augmentation  Stanton  number,  given  by: 


CH  u * 0.098  X°*3l7sin20 

G * local  body  angle 

X = PpVp(1+G)/pODUoo 

ppVp  * particle  freestream  mass  flux 
G * mass  loss  ratio 


This  correlation  is  applied  in  parallel  with  the  crater  roughness  heating  modeling, 
for  the  laminar  and  transitional  regions  of  the  nosetip.  That  is,  the  larger  of 
the  two,  usually  the  erosion  augmentation,  is  used  in  the  ablation  calculation. 


Particle  Slowdown/Ablation 

As  the  hydrometeor  traverses  the  bow  shock  layer,  it  is  exposed  to  relatively 
hign  temperature,  dense  air.  Under  these  conditions  water,  ice,  and  snow  particles 
slowdown  and  lose  mass  to  the  air.  The  model  developed  originally  in  Reference  18 
and  improved  slightly  i Reference  1 accounts  for  the  particle  slowdown  and  mass 
loss  but  does  not  consider  particle  deflection  in  the  shock  layer. 


Particle  Stripping  and/or  Breakup 

As  liquid  rain  drops  traverse  the  shock  layer,  data  indicate  that  mecha- 
nisms exist  lo  cause  surface  layers  of  the  drop  to  be  stipped  off,  thus  reducing 
the  effective  oarticle  mass.  Furthermore  at  certain  conditions  the  drops  can 
breakup  catastrophically.  Modeling  described  in  Reference  19  was  incorporated  in 
the  EROS  code  to  study  these  phenomena,  however,  the  residence  time  comes  from  the 
slowdown/ablation  calculation.  In  the  EROS  code  stripping  breakup  model.  Bond  number. 
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Bo*.  Is  assumed  to  be  In  the  range  where  stripping  and/or  breakup  will  occur 
(Bo  > 10).  As  such,  a critical  normalized  time,  T*  Is  defined  from  the  Weber 

t ^ 

number,  W , as  follows: 

Tc  * 45  W"°-JS 

For  T > T , breakup  occurs  and  no  mass  hits  the  surface.  For  T < T , no  breakup  occurs 
but  stripping  may  be  Important.  In  the  stripping  model,  the  mass  of  the  raindrop 
Is  given  by: 

So  * 7 (]  + cos  n rs) 


where 


M - impact  raindrop  mass 
Mq  = initial  raindrop  mass 


Note  that  for  T > 3.5  all  mass  is  stripped  off  and  the  effective  impact  mass  is  zero. 


For  the  one  sensitivity  solution  which  exercised  this  model,  the  procedure 
was  applied  to  all  hydrometeor  types  (water,  snow  and  Ice).  Furthermore,  only  the 
median  particle  size  was  considered.  The  effects  of  particle  size  distribution  in 
either  the  slowdown/ablation  or  the  stripping/breakup  calculations  were  not  addressed 
during  these  code  checkout  calculations. 


1 

Bo  = g ppFQ2/o 


Where  g = 


o = 


particle  acceleration 
particle  density 

particle  initial  radius 

particle  surface  tension 


t = tu2/oo(p2/Ppy/* 

Where  t = resident  k.ime  in  shock  layer 

u2  - gas  velocity  behind  shock,  relative  to  partic'e 
P2  = gas  density  behind  shock 
Dq  = initial  particle  diameter 

W = 2 P2u22r0/a 
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4.3  RESULTS  OF  EROSION  EFFECTS  SOLUTIONS 


Results  from  the  solution  matrix  (Table  4-2)  are  compared  In  this  section  to 
evaluate  the  sensitivity  of  nosetlp  response  to  erosion  modelling  variations;  these 
evaluations  are  presented  In  the  form  of  stagnation  point  recession  comparisons.  In 
addition,  flight  data  are  presented  as  ar  Indication  of  the  overall  modeling  accura- 
cy. The  ablation  modeling  for  graphitic  materials,  summarized  In  Section  3,  was 
used  In  the  erosion  effects  computations.  The  erosion  modeling  Is  described  In 
Reference  1 and  Section  4.2.  The  effects  of  respective  modeling  features  are  dis- 
cussed In  the  the  following  paragraphs. 

Baseline  Calculations 

Baseline  calculations  for  the  two  SAMS  flights  reveal  that  the  coupling  be- 
tween erosion  and  ablation  Is  significant.  Figure  4-4  shows  the  relative  contribu- 
tions to  the  total  recession.  Using  the  baseline  modeling,  only  50  to  60  percent  of 
the  recession  Is  directly  due  to  erosive  mass  loss.  The  ablation  component  (40  to 
50  percent  of  the  total)  results  primarily  from  heat  transfer  Increases  associated 
with  the  hydrometeor  encounter.  Clear  air  heat  transfer  and  ablation  produces  very 
little  nosetlp  recess lor  in  these  Terrier-Recruit  flight  tests  (less  than  0.1  In.). 
The  comparison  with  data  shots  that  the  baseline  modeling  gives  generally  too  much 
recession,  35  percent  high  fcr  SAMS  6,  16  percent  high  for  SAMS  7.  This  level  of 
differences  and  the  differences  between  the  computed  and  measured  recession  histor- 
ies are  well  within  cloud  profile  uncertainties.  The  solutions  described  below  pro- 
vide Insight  Into  erosion  reluted  modeling  which  significantly  affect  the  calculated 
recession  and  which  may  also  explain  the  differences  between  prediction  and  data. 

Minimum  Erosion  Mass  Loss  Correlation 

The  effects  of  uncertainties  In  erosion  mass  loss  data  were  assessed  by  using 
an  alternate  mass  loss  correlation.  Since  the  baseline  solutions  were  above  the 
flight  data,  a lower  bound  interpretation  of  the  mass  loss  data  given  In  Reference  1 
was  used.  As  described  In  Section  4.2, 


— = 4.44  V"0, 2 
horn 


where  V Is  the  velocity  of  the  Impacting  mass.  The  solutions  for  the  two  SAMS 
flights  are  compared  to  the  respective  baseline  computations  in  Figure  4-5.  The  re- 
sults are  a 20  percent  reduction  in  stagnation  point  recession  for  SAMS  6 and  22 
percent  reduction  for  SAMS  7. 
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(b)  S/VMS  7 


I Figure  4-4.  Baseline  eros lon/ablation  calculations. 

i 

I 
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(b)  SAMS  7 


Figure  4-5.  Erosion/ablation  calculations 

using  lower  mass  loss  correlation. 


For  the  velocity  conditions  typical  of  the  Terrier-Recruit  trajectory  (6  to 
8.6  kft/sec),  the  lower  "G"-1aw  gives  between  23  and  27  percent  reduction  In  mass 
loss  ratio  (G).  Hence,  the  erosion  component  of  recession  Is  reduced  directly  by 
these  amounts.  The  ablation  component  is  also  reduced  because  both  crater  roughness 
and  erosion  heating  augmentation  are  dependent  on  the  mass  loss  ratio.  The  overall 
20  to  22  percent  reduction  in  computed  stagnation  point  recession  consists  of: 

• Up  to  14  percent  from  reduced  erosion  mass  loss 

• About  7 percent  from  reduced  erosion  effects  on  heat  transfer  and  abla- 
tion 

These  reduction  do  provide  better  agreement  with  the  final  recession  measurements. 
However,  recession  histories  do  not  agree  as  well. 

Erosion  Effects  on  Heat  and  Mass  Transfer 

The  baseline  erosion  modeling  accounts  for  the  effects  of  erosion  augmenta- 
tion and  crater  roughness  on  surface  heat  and  mass  transfer.  In  the  ER05  code,  the 
modeling  of  these  phenomena  compete  with  the  clear  air  models  to  determine  the  heat 
transfer  boundary  condition.  The  relative  importance  of  cratering  effects  and  ero- 
sion augmentation  during  the  two  SAMS  flights  was  assessed  through  a series  of  EROS 
code  calculations  (Solutions  4 through  7 of  Table  4-2).  The  following  discussions 
compare  the  results. 

• Solution  4 shows  the  effect  of  complete  crater  healing  (i.i:.,  no  crater 
formation)  on  the  comDiited  recession 

c Solution  5 was  run  to  compute  the  amount  of  crater  roughness  reduction 
given  by  the  ablation  healing  model  derived  in  Reference  1 

• Solution  6 indicates  the  effect  of  eliminating  the  erosion  augmentation 
model . 

• Solution  7 allows  only  clear  air  heat  transfer  modeling  by  eliminating 
both  crater  effects  on  heat  transfer  and  erosion  augmentation.  Thus,  the 
solution  could  be  considered  uncoupled  (no  erosion/ablation  coupling). 

The  Solution  5 calculations  indicate  that  ablation  healing  does  not  signifi- 
cantly reduce  crater  roughness  for  the  ablation  rates  and  obscuration  times  typical 
of  graphitic  nosetips  in  a Terrier-Recruit  flight.  The  maximum  reduction  in  rough- 
ness was  roughly  10  percent.  Since  roughness  heating  factors  are  correlated  in 
terms  of  the  logarithm  of  roughness  (Reference  1),  the  10  percent  reduction  has  a 
negligible  effect  on  the  roughwall  heat  transfer.  Computed  recessions  for  SAMS  6 
and  SAMS  7 agreed  closely  with  the  baseline  results. 
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When  crater  roughness  Is  eliminated  entirely,  as  in  Solution  4,  clear  air 
scallop  roughness  dominates  the  turbulent  heating  computation,  and  erosion  augmenta- 
tion dominates  the  laminar  and  transitional  distributions.  As  demonstrated  in  Fi- 
gure 4-6,  however,  elimination  of  crater  roughness  has  a minimal  effect  on  the  com- 
puted recession  history.  Slight  differences  late  In  the  flights  are  associated  with 
shape  change  effects.  Similar  results  are  achieved  when  crater  roughness  is  allowed 
but  augmentation  is  eliminated  (Solution  6).  Comparative  recession  histories  are 
shown  in  Figure  4-7.  In  this  case,  craters  promote  transition  very  close  to  the 
stagnation  region  and  dominate  the  roughness  effects  on  laminar,  transitional,  and 
turbulent  heat  and  mass  transfer.  Since  the  crater  and  erosion  augmentation  effects 
are  modeled  to  act  in  parallel,  it  is  concluded  that  they  have  a nearly  equal  magni- 
tude for  these  flight  conditions. 

Figure  4-8  demonstrate*;  that  elimination  of  both  cratering  and  erosion  augmen- 
tation (Solution  7 of  the  sensitivity  matrix)  substantially  reduces  the  overall  com- 
puted recession.  The  recession  in  this  case  is  roughly  comparable  to  the  erosion 
plus  clear  air  ablation  line  shown  in  Figure  4-4.  Erosion/ablation  coupling  is, 
therefore,  a result  of  either  surface  cratering  or  erosion  augmentation. 

Figure  4-9  compares  the  stagnation  point  heat  transfer  coefficient  histories 
for  four  SAMS  7 calculations.  The  curves  demonstrate  the  significant  heat  transfer 
increase  when  either  erosion  augmentation  or  crater  roughness  are  included  in  the 
solution.  Furthermore,  since  the  no  crater  solution  is  identical  to  the  baseline 
case,  it  may  be  concluded  that,  for  the  Terrier-Recruit  environment,  the  erosion 
augmentation  controls  the  response  when  both  effects  are  being  considered. 

Particle/Shock  Layer  Interactions 

Two  particle/shock  layer  interaction  models  are  included  in  the  EROS  computer 
code.  Solutions  8 and  9 of  the  sensitivity  matrix  indicate  the  effects  of  these 
models  on  the  SAMS  6 and  SAMS  7 computations.  The  particle  slowdown/ablation  model 
significantly  reduces  stagnation  point  recession,  as  shown  in  Figure  4-10.  The  pri- 
mary cause  is  the  ablation  mass  loss  of  the  particle.  Very  little  slowdown  occurred. 
Indeed,  for  the  aft  cone  region  of  the  nosetip,  the  path  length  (shock  to  body)  is 
sufficient  to  allow  complete  particle  demise. 

Calculations  were  also  pe  formed  using  the  particle  stripping  and  breakup 
model.  This  nmdel  applies  only  to  raindrop  encounters.  Although  the  two  flights 
were  conducted  in  weather  containing  predominantly  ice  particles,  the  modeling  was 
exercised  for  the  full  trajectory  range.  The  median  particle  size  as  a function  of 
altitude  was  used  in  the  calculation.  In  all  cases,  the  modeling  had  a negligible 
effect  on  the  particle  mass.  Thus,  the  solutions  were  identical  to  the  baseline 
results.  This  is  reasonable  considering  the  particle  size  assumptions. 
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(b)  SAMS  7 


Figure  4-6.  Erosion/ablation  calculations 
allowing  no  crater  roughness. 
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Figure  4-8.  Erosion/ablation  calculations 
allowing  no  crater  roughness 
and  no  erosion  augmentation. 
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Figure  4-9.  Comparison  ,f  stagnation  point  heat  transfer 
coefficients  versus  altitude  for  SAMS  7. 
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(b)  SAMS  7 


Figure  4-10.  Erosion/ablation  calculations  allowing 
particle  slowdown/ablation. 


The  significant  results  of  the  weather  effects  sensitivity  calculations  for 
the  Terrier-Recruit  flight  conditions  are  as  follows: 


• A significant  portion  of  the  mass  loss  is  associated  station  due  to  the 
coupled  effects  of  hydrometeor  encounters  on  the  surface  transfer. 

• The  coupling  between  erosion  and  ablation  Is  accomplished  through  either 
crater  roughness  or  erosion  augmentation  modeling. 

• For  the  Terrier- Recruit  flight  conditions,  crater  roughnes:  and  erosion 
augmentation  effects  are  comparable  although  augmentation  generally  dom- 
inates slightly. 

• Ablation  healing  of  craters  reduces  crater  roughness  less  than  10  percent 
for  the  SAKS  6 and  SAMS  7 trajectory  and  weather  conditions;  this  reduc- 
tion makes  a negligible  change  in  surface  heat  transfer. 

• For  the  SAMS  6 and  SAMS  7 flight  cases,  the  particle  slowdown  and  abla- 
tion model  significantly  affects  predicted  recession  by  reducing  the  mass 
of  the  impacting  particles.*  The  particle  slowdown  was  small. 

• The  water  drop  stripping  and  breakage  model  had  little  effect  on  the  pre- 
dicted nosetip  response;  presumably  because  only  a single,  relatively 
large,  particle  size  was  assumed. 

• Because  both  erosion  and  ablation  at  Terrier-Recruit  flight  are  closely 
dependent  upon  the  erosion  mass  loss  ratio,  recession  is  directly  related 
to  changes  in  the  mass  loss  correlation. 

The  recovered  nosetips  from  SAMS  tests  indicate  both  the  overall  stagnation 
point  recession  and  the  nosetip  shape,  including  the  cone  surface  recession.  Com- 
parisons between  computed  and  measured  shape  profiles  show  the  adequacy  of  the  shape 
modeling  numerics.  Figure  4-11  presents  scaled  comparisons  of  shape  profiles  for 
SAMS  6.  The  measured  shape  is  shown  by  the  short  dashed  lines,  solid  lines  indicate 
the  two  solutions  from  the  sensitivity  matrix  which  best  matched  the  overall  reces- 
sion data  (Solutions  3 and  8 of  Table  4-2).  The  computed  side  surface  recession  is 
much  greater  than  observed,  even  for  the  slowdown/ablation  case  where  cone  region 
particles  demise  before  impact.  It  is  probable  that  transient  heat  conduction  ef- 
fects, not  modeled  in  the  sensitivity  solutions,  would  reduce  side  wall  recession 
substantially.  The  comparisons  of  nose  region  shapes  are  quite  favorable,  in  that 
both  the  measured  and  computed  were  quite  blunt. 


★ 

The  particle  slowdown/ablation  model  is  derived  tor  steady  state  particle  vaporiza- 
tion. Later  calculations  reported  in  PANT  Interim  Volume  XIX  showed  that  droolets 
do  not  vaporize  at  SAMS  conditions  (Mach  < 10.). 


L 
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Figure  4-11.  Relative  final  shape  comparisons  for  SAMS  6. 
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SECTION  5 


ALTERNATE-MATERIAL  EROS  CODE  CALCULATIONS 

Two  additional  SAMS,  Terrier-Recruit,  flight  cases  were  analyzed  using  the 
EROS  computer  code.  There  were  Sandia  Tests  Number  R487411  and  R487406.  The  objec- 
tive of  these  analyses  is  to  check  out  the  alternate  material  modeling  capability  of 
the  EROS  code  on  a charring  ablator  and  an  alternate  graphitic  material.  These 
tests  were  conducted  from  Wallops  Station,  Wallops  Island,  Virginia  on  8 April  1974. 
Trajectory  and  weather  data  were  obtained  from  Reference  16,  and  the  key  flight  test 
parameters  are  summarized  in  Table  5-1.  Pertinent  trajectory  quantities  are  shown 
in  Figure  5-1;  the  weather  parameters  (liquid  water  content  and  median  particle  dia- 
meter) used  in  the  computer  analyses  are  displayed  in  Figure  5-2.  The  nosetip  con- 
figuration is  unchanged  from  the  SAMS  6 and  SAMS  7 flight  (refer  to  Figure  4-2). 

TABLE  5-1.  SAMS  FLIGHT  TEST  PARAMETERS 


Parameter 

R487411 

R4 87406 

Peak  Velocity  ft/sec 

8600 

8600 

Peak  Stag.  Press.,  atm 

58 

58 

Nose  Radius,  inches 

.625 

.625 

Core  Angle,  degrees 

9 

9 

WSI,  km2-gm/m3 

2.1 

2.1 

Nosetip  Material 

Carbon/Carbon 

Carbon  Phenolic 

The  erosion  mass  loss  modeling  for  the  carbon  phenolic  material  (Flight 
R487406  is  summarized  in  Reference  20.  The  mass  loss  correlation  for  Flight  R487411 
taken  to  be  the  one  recommended  for  use  with  h III -A  carbon/carbon  in  Reference  1. 
It  has  the  form: 


G oc  Visin'  ) 

The  ablation  response  modeling  of  carbon  phenolic  is  reviewed  in  Reference  20.  The 
modeling  of  the  carbon/carbon  was  assumed  identical  to  that  used  for  the  SAMS  6 and 
SAMS  7 solutions  (Section  4). 
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The  analysis  results  are  presented  as  shape  profile  histories  and  stagnation 
point  recession  histories  In  Ftgures  5-3  and  5-4. 

The  conclusions  Is  that  the  EROS  code  functioned  properly  for  the  two  alter- 
nate material  flight  test  cases.  No  data  were  available  at  this  writing  for  compari- 
son with  the  EROS  code  solutions. 
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(b)  Recession  history 


Figure  5-3.  Nosetip  response  prediction  for 

SAMS  Flight  R487411,  carbon/carbon. 
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SECTION  6 


CONCLUSIONS  AND  RECOMMENDATIONS 

Conclusions  are  sumnarized  in  Section  6.1  and  recommendations  are  listed  in 
Section  6.2. 

6.1  CONCLUSIONS 

The  EROS  code  was  exercised  to  determine  code  reliability  and  computational 
consistency  for  wind  tunnel  cest  environments,  clear  air  fliqht  environments,  and 
weather  flight  env.  >nments.  Solutions  were  performed  exe  ing  various  modeling 
perturbations,  and  from  these  solutions,  the  sensitivity  of  the  EROS  code  predic- 
tions to  key  modeling  variations  and/or  uncertainties  was  demonstrated.  The  primary 
conclusions  from  the  code  calculations  may  be  summarized  as  follows: 

Comparisons  to  Mind  Tunnel  Data 

1.  Slight  improvements  in  the  calculations  of  bow  shock  and  turbulent  heat 
transfer  were  achieved. 

2.  Shape  change  solutions  compared  favorably  with  low  temperature  ablation 
data,  although  for  one  case  overall  recession  was  inaccurate  because  of 
deficiencies  in  scallop  roughness  modeling. 

Clear  Air  Shape  Change  Sensitivity  Studies 

3.  The  total  recession  of  an  ICDM  nosetip  at  typical  conditions  is  over- 
whelmingly dependent  upon  the  transition  altitude. 

4.  Based  on  the  PANT  roughwall  ’•‘insition  criterion,  transition  altitude  is 
extremely  sensitive  to  the  vuiue  of  the  surface  roughness  which  exists 
during  laminar  ablation;  greatest  sensitivity  occurs  for  roughness  values 
between  0.2  and  0.4  mils  for  typical  reentry  situations. 

5.  If  the  PANT  roughwall  heating  model  is  correct  and  applicable  to  ablation- 
induced  scalloped  surfaces,  then  a factor  of  two  de;rease  in  scallop  di- 
mension produces  a relatively  small  (12  to  18  percent)  decrease  in  nose- 
tip  recession  rate. 
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6.  The  transitional  heat  transfer  modeling  affects  the  flight  time  required 
for  nosetip  sharpening. 

7.  Two-dimensional  heat  conduction  is  important  but  not  critical  to  the 
evaluation  of  shape  and  recession  trends  for  graphitic  materials.  The 
two-dimensional  heat  conduction  towards  the  centerline  generally  in- 
creases stagnation  point  recession  and  produces  a slightly  larger  nose 
radius  shape  than  is  predicted  ass-ming  steady-state,  one-dimensional 
heat  conduction. 

Coupled  Erosion/Ablation  Sensitivity  Calculations  for  Graphite 

8.  For  the  Terrier-Recruit  flights  analyzed,  the  coupled  effect  of  erosior, 
on  ablation  mass  loss  is  very  significant.  Computed  ablation  including 
weather  effects  is  40  to  50  percent  of  total  mass  loss.  Coupled  effects 
would  also  be  significant  for  light  or  mild  weather  encounters  at  ICBM 
conditions . 

9.  The  hydrometeer  encounter  increases  surface  heat  transfer  and  ablation 
through  either  erosion  augmentation  or  crater  roughness  formation.  For 
the  Terrier-Recruit  conditions,  erosion  augmentation  was  slightly  greater 
than  the  crater  effects. 

10.  Uncertainties  in  mass  loss  ratio  correlations  have  a significant  impact 
on  the  nosetip  recession  calculation.  However,  mass  loss  uncertainties 
at  Terrier-Rec  uit  conditions  are  probably  within  weather  definition  un- 
certainties. 

11.  For  the  nominal  erosion  mass  loss  correlation,  modeling  of  particle  slow- 
down and  ablation  in  the  shock  layer  significantly  affected  the  computed 
response. 

12.  The  crater  healing  and  particle  stripping/breakup  models  had  negligible 
effect  on  the  Terrier-Recruit  calculations. 

Alternate  Material  Erosion/ftblation  Calculations 

13.  Two  coupled  erosion/ablation  shape  change  solutions  were  performed  to 
check  out  the  generalized  material  modeling  capability  of  the  EROS  code. 
These  solutions  were  successfully  completed. 
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6.2  RECOMMENDATIONS 


Recommendations  apparent  from  the  EROS  code  checkout  calculations  are  as 
.follows: 

1.  Since  clear  air  nosetip  response  is  critically  dependent  upon  the  intrin- 
sic surface  roughness  during  laminar  ablation,  this  material  property, 
its  uncertainty  and  its  variability  within  a particular  billet  should  be 
thoroughly  characterized  prior  to  qualification  of  the  material. 

2.  Additional  effort  is  required  to  understand  the  effects  of  roughness  var- 
iabilities on  the  nosetip  response. 

3.  Applicability  of  the  PANT,  roughwall,  turbulent  heating  model  to  scal- 
loped surfaces  should  be  assessed  in  relation  to  other  techniques  cur- 
rently in  use. 

4.  The  sensitivity  of  EROS  code  predictions  to  time  step  size  and  surface 
coordinate  definition  procedures  should  be  evaluated. 

5.  The  EROS  code  heat  conduction  procedure  requires  generalization  to  allow 
multiple  in-depth  back-up  materials  and  improvement  of  stagnation  point 
recession  calculation. 

6.  Physical  justifications  for  the  erosion  mass  loss  correlations  and  the 
erosion  augmentation  formulations  in  current  use  are  required  before  ex- 
trapolation to  I CBM  conditions  can  be  done  with  confidence. 

7.  Particle/shock  layer  interactions  should  be  included  in  coupled  effects 
modeling  for  Mach  > 10.* 


See  PANT  Interim  Volume  XIX. 
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APPENDIX 

SHAPE  PROFILE  HISTORIES  FROM  CLEAR  AIR  PARAMETRIC  SOLUTIONS 

Table  A-l  identifies  the  fifteen  shape  response  solutions  completed  as  a 
part  of  the  EROS  code  check  out  calculations.  The  shape  profile  histories  for 
these  solutions  are  given  in  Figures  A-l  through  A-15. 
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Solutions  9 and  12  performed  using  both  transient  and  steady  state  conduction  modeling. 


Figure  A-2.  Solution  2,  case  1 smeller  microroughness. 


•e  A-3.  Solution  3,  case 
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Figure  A-10.  Solution  9,  case  2 alternate  transition  criterion  (transient). 
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Figure  A-12.  Case  2,  alternate  transition  criterion 
plus  lower  macrorougfwess. 


